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Travelling To and Through Space 


Space Launch System (SLS) - 
America’s Heavy-lift Rocket 


Provides initial lift capacity of 70 metric tons (t), 
evolving to 130 t 


Carries the Orion Multi-Purpose Crew Vehicle 
(MPCV) and significant science payloads 


Supports national and international missions 
beyond Earth’s orbit, such as near-Earth asteroids 
and Mars 
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How turbomachinery is used in Rocket Engines 



are stored in Fuel tanks at a few atmospheres. 

• Turbines, driven by hot gas created by mini- 
combustors, tied with shaft to pump, which 
sucks in propellants and increases their 
pressures to several hundred atm. 

• High pressure propellants sent to Combustion 



chamber, which ignites mixture with injectors 

• Very hot gas directed to converging/diverging 
Nozzle to increase flow to very high velocity for 
thrust. 



Motivation: Avoid High Cycle Fatigue Cracking in Turbomachinery 


• Cracks found during ground-test program stop engine development 

- If cracks propagates, it could liberate a piece, which at very high rotational speeds 
could be catastrophic (i.e., engine will explode). 



• In J2-X Rocket Engine program, became apparent that 
turbine blade external damper (needed to show 


deterministic design good) behind schedule. 


Identified probabilistic analysis as method to quantify risk 


during individual tests in series. 





Prediction of Probability of Failure of Turbine Blades during 
Testing - Motivation 

• Standard blade forced-response analysis process recognizes uncertainty in 
material properties and in prediction of natural frequencies. 

• For J-2X clear that other non-deterministic variables (damping, mistuning) also 
important. 

• Needed to properly assess risk of blade failure using actual non-deterministic 
nature of these rv's rather than using deterministic design values. 

• Substantial research and application in literature of probabilistic methods to 
turbomachinery issues 

- Much of effort ("top down") calculates reliability by comparison to measured 
reliability of sub-systems on similar engines - Packard, '02. 

- Crack growth characterization in probabilistic FEA ("bottoms-up) - Petrov, '08. 

• OBJECTIVE -calculate probability of failure using closed-form finite life solutions 
in terms of these 4 non-deterministic variables and peak FEA-derived stress 
state. 

• Answer 1) What is P f during a specific test series? 

2) If previous analysis showed low safety factors, why didn't it fail? 


t 


• Variation of Natural Frequency, fn 
typically accounted for using rule- 
of-thumb +/-10% in frequency 
response analysis. 


• Here, data from previous engine 
programs show distribution is 
somewhat Gaussian with a 3a 
variation of +/- 5% (p=1.67%). 



MocJe 2r. 8. Slandaifi beuauon =i74.i i 

IE. ODD -1 — 





Input Variables & Assumptions 


• In gas-generator rocket engine cycle flow rate, turbine efficiency determines rotor 
rotational rate, so Operating Speed is random variable. 

• For the engine balance used here, resulting operating speed distribution is 


Speed ~ N(|j=30,635 rpm, g= 307.7 rpm) 
^ +/- 3a range is 6% 


• Exception is in "Powerpack" testing, where turbopumps are isolated and 
rotational speed is controlled. 


Input RV’s - Mistuning Background 

• Imperfectly cyclically-symmetric (mistuned) bladed-disks exhibit warping of nodal- 
diameter modes and amplification of peak response compared to perfect cyclic- 
symmetric (tuned) disk. 

• Effects of mistuning non-deterministic since every build will be different. 

• J-2X is one of first rocket turbopumps developed since practical methods developed 
to predict statistics of mistuning amplification value m. 

• Analysis performed ("SNM method") to develop statistics of m for 3 of J-2X 
problematic modes. 





Statistics of Amplification due to Mistuning 


• For 69 blade-disks, stats developed for entire-blade population (690,000) 
and max-responding blade per bladed-disk population (10,000) for 3 
different modes. 

• Debate -> concensus: for probabilistic analysis, use mean value of 0.9 with 
Lognormal fit. 



m 



Input RV’s - Damping 


• Damping is critical parameter for forced response prediction, so "whirligig" test 
program used to obtain data. 

• Whirligig was mechanically-driven rotor with similar bladed-disk (J2-S) with 
similar dampers excited by pressurized orifice plate to simulate blade excitation. 

• Key assumption is that this reflects true configuration. 


• SDOF Curve fit 
technique applied to 
selected top-responding 
blades to derive 
damping from response. 




Damping Results from Whirligig 


Data shows wide-variation in 
damping, but reasonable 
population (15-20 acceptable 
samples) for characterization 
of mean and type. 

Lognormal distribution fits 
obtained for each mode. 
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Probabilistic Analysis 


• First, determine Stress state (Sa, Sm) of problem location from finite 
element frequency response analysis at resonance (w. ^=.0025). 

• Then, for a sample taken from distributions of all random variables 
(ie, Monte Carlo analysis), calculate Equivalent Alternating Stress A eq : 
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Nominal HCF cycle count data ("s-n curves") -> 

N _ JQ(- 9 - 2461xio ^io(4 ? )+20.672) 


Finally, failure using "finite life approach" is 
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Probability of Failure using Damage Fraction 


N 


accum 


is excitation frequency * time at that frequency. 


- Speeds in test series recorded in 80 rpm wide bins, calculate incremental damage 
fraction within each /'th bin. 
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Technique Verification 


• To "verify" technique, p f was calculated for tests that had already taken place, 
assuming both that the speeds are "post-priori" known and "a-priori" 
unknown. 

• Deterministic analysis indicated Safety Factor < 1 for mode 14 in ND 5 family 

- fn ~ N(36851 hz, 615 hz), £ ~ LN(0.304%, .097%). 

• Results for these technique verifications were reasonable 

- for a single hot-fire test, p f only 1% (specifically because of a low 
probability of resonance) , so fact that blade did not crack should be 
expected. 


J-2X Powerpack II A-priori Analysis 


• Speed mean controllable, enabling engineering team to make assumptions of 
1200s total run time, 4 dwells of 100s, 20 dwells of 30s, ramping from 26902 rpm- 
31200 rpm at continuous rate of 20 rpm/s (during ramps). 

P f = 33.8% 

• Explanation for results is extensive overlap of fn distribution over operating range 
(ensuring resonance) , and lower damping of problematic mode. 



Operating Range (rpm) 


• Test Results - dampers not put in, extra precautions taken, blade did not fail 
"Statisticians are never wrong , they are only unlucky" 





J-2X Engine 10001 A-Priori Analysis 


• Equally important to assess P f for first full-scale engine test to 
determine if external blade dampers required. 

• In this test speed will resolve to a single value within distribution 
following Speed ~ N(30635 rpm, 307 rpm). 

• Time of operation given as 550 s. 

• Single dwell formulation relatively simple, enables large (100,000) 
sample MC run. 

• P f = 1.06%, very low because of low probability of resonance itself, 
which was independently calculated (using only rv's speed and 
natural frequency) to be P resonance = 3.1%. 


Sources of Error and Conclusions 


• Error: 

- Some non-deterministic input variables assumed to be deterministic. 

- Mistuning and Damping assumed to be independent and they probably are not; 
unknown effect on results. 

- Response away from resonance approximated by SDOF curve fit. 

• Framework procedure established for quantifying risk of turbine blade failure due 
to resonance. 

• Probabilistic analysis enable first-time use of statistical distributions of most of 
random variables, including Natural Frequency, Operational Speed, Mistuning, 
and Damping. 

• Results very useful for project decision-making during development phase. 

• Framework also applied to a number of other J2X turbopump dynamics issues. 

- Used to determine appropriate deterministic value of damping to use for design for 
specific reliability goals. 

- Design of test series to put equivalent damage on pump inducer blade as it would 
experience if it were at resonance (worst case), given that the fn is actually non- 
deterministic. 





Combination of Random and Harmonic Loads in 

Structures - Introduction 


• Many structural components are in an environment with both random and 
harmonic loads. 


• Rocket Engines 

Turbopump-harmonic 


Combustion-random 


• Each type first 
calculated 
individually. 


• Results of analyses 
then combined for 
use by stress in both 
ultimate/yield analysis 
and HCF analysis. 


e.g., 1 lb Sine 
Amplitude load at 1 hz 


Frequency response analysis to 
generate harmonic load first 
calculated n_ I 


t (sec) 




Random Loads 



Mean Square ® 2 =Area under random response curve 
= 400 lb 2 

lcr = RMS = Vo 1 = 20 lb 


Loads extremely sensitive to probability level 
chosen (or actually obtained) because of flatness of 
Cumulative Distribution Function at tails. 






Loads Combination Equations 


• Extensive, difficult research into reducing each load type individually; however, 
little thought into how load is combined. 

• Main goal of methods is to estimate an "equivalent 3a" design load; 

- 3a is traditionally used for pure random loading, i.e., the load that exceeds 
99.865% of the occurrences. 

• " Standard Method " used in SSME: 

design load = A sm + 2>a mn 


• "3*ssMS" Method: 


design load = 3 J(<r sm ) 2 +(a rando J 


• Both techniques exceed 99.865% by definition, not tied to a specific probability. 

• " Peak " method proposed by Steinberg, adopted initially by engine contractors. 


design load = ^J(A sin ) 2 +(3cr mni!om ) 2 




Typical MC-1 Engine Load Set 
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SB-6 
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Loads Combination using PDF’s 


• Harmonic signal can be defined as stationary random process when combined 
with an independent Gaussian process since phase relationship with random 
signal is random. 

• Define harmonic signal as 

•Vsin = AsVd((Ot + (p) 

• Then PDF of sine distribution is 




“Exact” Solution Now Easily Obtained 


Create and Integrate Joint PDF of Normal and Sine Distributions to obtain 
CDF of design load z: 
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Mathematico® can perform not only integration, but also inverse: 


- Given a load (e.g. calculated using "standard" method) calculate exact 
reliability level. 

- Given a desired reliability level, solve for corresponding load. 


Developed Excel Macro: 

- easily integrates into existing loads calculation spreadsheets 

- Accesses Mathematica® to perform inverse-integration to obtain 
design load corresponding to 99.865% reliability 

- returns value seamlessly into spreadsheet. 


Loads Combination using Monte Carlo 

• Gaussian random vector using random analysis results (a ran ) first simulated: 

{r}~N(0.0,o random ) 

• Independent sine vector generated using harmonic analysis results (Aj): 

- Create uniform distribution {x}j~ U(0,1) 

- Generate sine distribution {y}j = Aj sin(2 n {x}j ) 

• Vectors of same length added to form total response: 

{z} = {r} + {y} 

• CDF calculated for {z}, 99.865% (or any other desired level) selected. 

• Excel Macro created to perform Monte Carlo Simulation to obtain design 
load corresponding to 99.865% reliability (within Excel). 

- Less than a minute for 400,000 samples. 




Example 







Evaluation and Comparison of Methods 
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Sine 

1 a 

random 

Integration 
method for 
99.865% 
(baseline) 

MC 

99.860 

(400,000 

samples) 

standard 
method - 
A + 3 ct 

% over- 
shoot 
from 
baseline 

CDF value 
from 

integration 

3*ssMS 

% over- 
shoot 

CDF value 
from 

integration 

"peak" 
method - 
ss(3sig,A) 

% over- 
shoot 

CDF value 
from 

integration 

10 

5 

22.034 

22.031 

25 

13.5% 

99.970% 

25.981 

17.9% 

99.990% 

18.028 

-18.2% 

93.930% 

5 

5 

17.668 

17.653 

20 

13.2% 

99.957% 

18.371 

4.0% 

99.912% 

15.811 

-10.5% 

94.896% 

5 

20 

60.915 

60.888 

65 

6.7% 

99.919% 

60.930 

0.03% 

99.865% 

60.208 

-1.2% 

95.758% 

26 

4 

34.760 

34.772 

38 

9.3% 

99.983% 

56.445 

62.4% 

100.000% 

28.636 

-17.6% 

94.291% 

97 

14.67 

129.081 

129.195 

141.01 

9.2% 

99.986% 

210.422 

63.0% 

100.000% 

106.517 

-17.5% 

94.316% 

50 

98.7 

313.047 

313.422 

346.1 

10.6% 

99.951% 

314.524 

0.5% 

99.871% 

300.292 

-4.1% 

95.534% 

64 

109.33 

352.240 

353.079 

391.99 

11.3% 

99.955% 

354.978 

0.8% 

99.875% 

334.176 

-5.1% 

95.443% 


• MC closely agrees with Integration method 

• Two generally accepted methods always above 99.965%. 

• “Peak” method underpredicts “3 a” value 





Curve Fit of Overshoot of 3*ssMS Method over CDF of 99.865% 
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• Similar equation derived for "Equiv. 
2a" (97.725%, research suggests more 
appropriate for HCF) 



Conclusions 


• Probability Values calculated, compared, & evaluated for several 
industry-proposed methods for combining random and harmonic 
loads. 

• Two new excel macros written to calculate combined load for any 
specific probability level. 

• Closed form Curve fits generated for widely used 3o and 2 a 
probability levels. 

• For design of lightweight aerospace components, obtaining 
accurate, reproducible, statistically meaningful answer critical. 




Accounting for Speed Variation (Dither) of Turbomachinery in 
Analysis - Introduction 


• Structural (S u!t & HCF) assessment critical for turbomachinery flow 
path components undergoing possible resonance. 

• Resonance generally avoided, but impossible for higher modes 
found with modern analysis, especially with wide speed ranges. 

- J2-X Fuel Pump turbine stator operates from 26Krpm-34Krpm; 69N forcing 
excites modes 10-18 between 30KHz-40Khz. 

• Criteria triggers forced response analysis at worst case resonant 
condition. 

• Finite life analysis, where actual fatigue damage during 
operational time is calculated, frequently used if endurance limit 
criteria violated. 
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Many Turbopumps “Dither” 


May be beneficial to incorporate fact that real turbopumps dither about a 
nominal mean speed, (separate from uncertainty in mean speed itself) 



J2-X Powerpack 
Adjusted Speed 
Trace 


During time speed is not exactly at natural frequency, damage accumulation is 
significantly reduced. 




Literature, Purpose 


• Initial studies of response of systems with time varying excitation frequency Q 
by Lewis- 1932, Cronin- 1965. 

• Lollack, 2002, defined reduction in peak response for monotically varying Q, 
useful for defining rate of sine-sweep tests. 

• Henson, 2008, studied harmonically varying Q. 

• For rocket engines, Q varies non-deterministically. Motivated previous work by 
authors (2010) that developed numerical approach for calculating response and 
general sensitivities. 

• Unacceptable HCF factor for J2-X stator resonant 30Khz mode prompted need 
for practical technique. 


• Purpose of this research 

- to develop practical design techniques that account for 
excitation frequency stochasticity in the fatigue life of 
turbomachinery components. 




Excitation Data 


' j 


• Taken from hot-fire testing of 
J2-X and SSME. 

• Q = engine speed 
(hz)* [forcing pressure 
distortions/Rev] (FPR). 

• Since purpose is to examine 
fatigue life at resonance, 
actual mean speed adjusted 
to natural frequency for 
analysis. 

• Histograms for two different 
engines show ~ Gaussian 
distribution of speed. 





Theoretical Basis, Numerical Transient Solution 


SDOF EoM 
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f(t) = Asin(m) 


Q is derivative of (|)(t), constant in classical vibration analysis. 

For specified time-varying Q, t 

m = | Cl(r)dT 
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Calculate A necessary to generate peak resonant value of a a)t previously 
obtained by FEA, 

A 
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CO 2 L, 


Now can solve for a a)t in EOM with using numerical Runge-Kutte procedure 
implemented in Matlab; agrees with Lollack's results for linearly varying Q. 


Finally, Calculate damage fraction ® using Miner's rule, ® = 
becomes 
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Peak Dynamic Stress (psi) 


Convergence of Time Step in Transient Solution 


• Applied deterministic speed variation from specific hot-fire test. 

• Time histories of Peak Dynamic Stress and Damage Fraction generated. 

• Convergence studies performed -> At=l/120f n . 




Time (s) 





Analytical Solution 




d (^sp ccd^ 

Hypothesis from previous work that if /» D — ~r — ; then closed-form 
(computationally fast) standard analytical equation for SDOF steady-state 


response would be accurate. 



• Validation by comparing 
response with numerical 
transient solution. 


Stress Response Zoom-In 





“Dither Life Ratio” for Specified Excitation History 


Calculation of damage performed considering dither for specific 10 sec. 
window. 

Damage calculation assuming constant resonant excitation -> 2.135 times 
more damage, call it "Dither Life Ratio". 




During design phase, actual speed time histories unknown, but statistics from 
similar engines known. 

Prompted development of Monte Carlo method using rapid analytical solution. 

Speed vector created using Normal 
distribution. 

Powerpack data -> std dev =38.6 hz 
(cov=0.129%). 

MC results linear because rate of 
change of frequency variation not 
correct (and very high), but damage 
accumulation is accurate on the 
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Sensitivity of DLR to speed COV and C, 


Accuracy of Monte Carlo technique with analytical solution allows comprehensive 
sensitivity study to key parameters 

Results: Larger for high COV for speed, since more time spent off- resonance. 

- Larger for small since peaks are sharper and time spent off-resonance will 
have less response. 


Blue Points - data 
Red Curves - fits 
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Conclusions 


• Numerical and Analytical methods developed to determine damage 
accumulation in specific engine components when speed variation 
included. 

• Dither Life Ratio shown to be well over factor of 2 for specific example. 

• Steady-State assumption shown to be accurate for most turbopump 
cases, allowing rapid calculation of DLR. 

• If hot-fire speed data unknown, Monte Carlo method developed that 
uses speed statistics for similar engines. 

• Application of techniques allow analyst to reduce both uncertainty and 
excess conservatism. 

• High values of DLR could allow previously unacceptable part to pass 
HCF criteria without redesign. 

• Given benefit and ease of implementation, recommend that any finite 
life turbomachine component analysis adopt these techniques. 


